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Abstract 

Using matrix-model methods we study three different N = 2 models: \J(N)xXJ(N) 
with matter in the bifundamental representation, U(iV ) with matter in the symmetric 
representation, and U(iV) with matter in the antisymmetric representation. We find 
that the (singular) cubic Seiberg-Witten curves (and associated Seiberg-Witten dif- 
ferentials) implied by the matrix models, although of a different form from the ones 
previously proposed using M-theory, can be transformed into the latter and are thus 
physically equivalent. We also calculate the one-instanton corrections to the gauge- 
coupling matrix using the perturbative expansion of the matrix model. For the U(iV) 
theories with symmetric or antisymmetric matter we use the modified matrix-model 
prescription for the gauge-coupling matrix discussed in ref. [1]. Moreover, in the ma- 
trix model for the U(iV) theory with antisymmetric matter, one is required to expand 
around a different vacuum than one would naively have anticipated. With these mod- 
ifications of the matrix-model prescription, the results of this paper are in complete 
agreement with those of Seiberg-Witten theory obtained using M-theory methods. 



1 Introduction 

The study of M = 2 gauge theories using the matrix-model techniques of Dijkgraaf and 
Vafa was initiated in [2, 3] (see [4] for earlier work) and was further developed and studied 
in [5, 6, 7, 8, 9, 10] (see also [11]). In [7] we showed that it is possible to derive all the building 
blocks of the Seiberg-Witten (SW) solution [12] (i.e. the curve and a preferred meromorphic 
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differential) purely within the matrix-model context. Since the curve is obtained from the 
large-M solution of the matrix model, one can obtain the SW curve in this manner only 
when an explicit large-M solution is available. However, as was stressed in [3, 5], even when 
the large-M solution is not available, one can still resort to perturbation theory to derive 
the prepotential order-by-order without knowledge of the curve or differential. A necessary 
ingredient for this procedure is the knowledge of the quantum order parameters Oj (periods 
of the SW curve); in [5] we proposed a perturbative definition of the periods. 

In this paper, we extend the matrix-model program to three Af = 2 gauge theories whose 
SW curves are non-hyperelliptic: U(iV) xU(iV) with a bifundamental hypermultiplet, U(iV) 
with a symmetric hypermultiplet, and U(iV) with an antisymmetric hypermultiplet. In sec. 2 
we review some previous results for these theories and reformulate them in a way that will 
facilitate comparison with our later discussion. 

The equivalence of the Af = 1 versions of the gauge theories considered in this paper and 
the corresponding matrix models was established, following the approach in [13], in refs. [1, 
14] (related earlier and subsequent developments can be found in [15, 16, 17, 18, 19, 20]). 
In principle these general results also show the validity of the matrix-model approach in the 
M = 2 limit. However, to obtain precise information about the Af = 2 gauge theories requires 
substantial additional work. Furthermore, there are several subtle issues, particularly for the 
theory with antisymmetric matter, that need to be addressed in order to recover the known 
results from the matrix models. 

In sec. 3 we compute the gauge-coupling matrices r^- using the perturbative expansion 
of the matrix model. For the XJ(N)xU(N) model, the gauge-coupling matrices are given by 
the second derivative of the free energy of the matrix model [21, 15, 13]. However, for the 
U(iV) theories with symmetric or antisymmetric hypermultiplets, certain crucial signs must 
be included among the terms of the second derivative of the free energy to obtain r^. The 
rationale for this was described in ref. [1], and is implemented in sec. 3.2.2. 

Furthermore, for the U(iV) theory with antisymmetric matter, we will show in sec. 3.2 
that, in order to obtain a gauge-coupling matrix and Af = 2 prepotential that agree with 
those computed using SW theory, one must perturbatively expand the matrix model around 
a vacuum different from the one that would have been naively anticipated. In subsequent sec- 
tions, we will see that this choice of vacuum is similarly required to reproduce the known SW 
curve and differential. The underlying reason for this choice of vacuum remains somewhat 
obscure to us. 

The matrix models associated with the gauge theories we consider can be solved in the 
large-M limit, giving rise to cubic algebraic curves [2, 23, 24, 25, 1] 

u 3 - r{z) u - s{z) = 0. (1.1) 

The functions r(z) and s(z) are determined by the matrix-model potential, up to two arbi- 
trary polynomials r± (z) and t\ (z) (which can be expressed in terms of the eigenvalues of the 
adjoint field(s) of the matrix model [24, 25, 1]). To fix the forms of r\(z) and ti(z), one must 
impose an additional criterion, namely, the extremization of the effective superpotential of 
the associated gauge theory. 

In sec. 4, we use matrix-model perturbation theory to provide a simple and efficient 
method for determining the polynomials r\(z) and ti(z), and thus the cubic curve, in an 
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expansion in the quantum scale A. We show that the cubic curves obtained for each of the 
theories considered can be transformed into the SW curves of those theories obtained using 
M-theory (and also from geometric engineering for the U(iV) xU(iV) theory). 

Also in sec. 4, we use extremization of W c g together with the saddle-point solution to 
derive a condition which implies, via Abel's theorem, the existence of a certain function on 
the matrix-model curve. We then show that such a function exists on the known (exact) 
M-theory curves. Assuming uniqueness, this demonstrates that extremization of W c g leads 
to a matrix-model curve that agrees exactly with the M-theory curve. 

The Seiberg-Witten differentials for the gauge theories studied in this paper can also be 
obtained within the matrix-model framework. In sec. 5, we compute these using matrix- 
model perturbation theory, obtaining agreement with the SW differentials known from M- 
theory. 

Appendix A contain a derivation of the SW curve and differential for the M = 2 \J(N) 
theory with fundamental hypermultiplets using methods developed in sec. 4 and 5 of this 
paper. Appendix B contains some technical details of the calculations of section 3. 



2 Cubic Seiberg-Witten curves from M-theory 

The Seiberg-Witten curves and differentials for the M = 2 gauge theories considered in this 
paper were previously obtained using M-theory methods, following the approach of ref. [26]. 
(The U(iV) xU(A r ) curve can also be obtained using geometric engineering [27].) 
The SW curves for the theories: 

(a) M = 2 U(iV)xU(iV) with an M = 2 bifundamental hypermultiplet [26, 28], 

(b) J\f = 2 U(iV) model with one symmetric (m) hypermultiplet [29], 

(c) J\f = 2 U(A r ) model with one antisymmetric (Q) hypermultiplet [29], 
obtained from M-theory considerations are given by 

y 3 + P(z) y 2 + A' N P(z) y + (A' 2 A'f = , (2.1a) 

y 3 + P(z) y 2 + A' N - 2 z 2 P(-z) y + A' 3N ~ 6 z 6 = , (2.1b) 



y 3 + 



„, , 3A ,N + 2 
P(z) + 



y 2 + 



A /JV+2 r , , 3A' N+ 
P(-z) + 



z 2 



z 2 



A / 3iV+6 

y + — — = o, (2.ic) 



z" 



where P(z) = Y\^ = i(z — and P(z) = Y\iLi(z — eQ- The Seiberg-Witten differential for each 
of the theories above, obtained from the M-theory setup, is given by [26, 30, 29] 

Asw — z — . (2.2) 

y 

The map y — > (A'A') N /y in the curve (2.1a) corresponds to exchanging the two factors 
of the gauge group, i.e. it can be undone by interchanging e[ <-> e[, A' <-> A', and thus 
leads to a physically equivalent curve. The two curves (2.1b), (2.1c) are invariant under 
the involutions z — > —z, y — > A' 2N ~ 4 z 4 /y, and z — > —z, y — > A' 2N+4 /z 4 y, respectively. 
The actual SW curves for these theories are the quotients of the curves (2.1b), (2.1c) by 
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the involution. This reflects the presence of the orientifold plane in the type IIA brane 
configurations that lift to the M-theory configurations leading to these curves. 

Using (2.1a)-(2.1c) and (2.2), the leading term in the instanton expansion of the pre- 
potential for each of the theories above was derived in refs. [31, 32, 33]. Recently a more 
efficient method has been developed [34, 35] based on earlier work [36]. (See also [37] for an- 
other approach.) In sec. 3, we will reproduce these results from a perturbative matrix-model 
calculation. 

For later comparison with matrix model results (sec. 4), we need to transform the above 
curves into another form, which is invariant under the maps discussed above. To do this we 
define 4 



u 



w = \z-V 



A / 2JV-4/3 _4/3 

y - 2A' N - 2 ^z^ - — , (2.3) 

V 

(2.4) 



y 



y 



where (5 — 0, 1, —1 for curves (2.1a), (2.1b), and (2.1c) respectively; note that u' is invariant 
under the maps discussed above. The variables u', w are related via 

Aw 2 = u' 2 + AA' N u', (2.5a) 
4z 2 w 2 = u' 2 + 4A >N - 2 z 2 u', (2.5b) 
W = z 2 u' 2 + 4A' N+2 u' . (2.5c) 

Using eq. (2.3) and eqs. (2.1a)-(2.1c), one may show that 



v!_-P(z) + 3A' N \ 

7 

y = -a^ t: j ^^i, i, (2.6b) 



P(z) + 3A' N - 2 z 2 



y - — hr^irl- (2 - 6c) 



Next, use the definition of v! to write w = z Py + ^z @u' + A' N 2f3 z 13 . Substitute eqs. (2.6a)- 
(2.6c) into this equation, square it, and use eqs. (2.5a)-(2.5c) to find 5 

iu' - P(z) + 3A' N ) v! («' - P(z) + 3A' N ) = A' N (P(z) - P(z)) 2 , (2.7a) 
(u' - P{z) + 3A' N - 2 z 2 ) u' {u 1 - P(-z) + 3A' N ~ 2 z 2 ) = A ,N - 2 z 2 (P(z) - P(-z)) 2 (2.7b) 

A/JV+2 

W - P(z)) u' (u' - P(-z)) = —^(P(z) - P{-z)) 2 . (2.7c) 

z z 

In sec. 4, we will compare the matrix-model curves to the SW curves written in this form. 



4 Hcnceforth, for simplicity, we will set A' = A' in the U(iV) xU(AT) model. 

5 The form of the curve for U(AT) with one antisymmetric hypermultiplet (2.7c) was first obtained in sec. 7 
of ref. [29], where the connection to the Atiyah-Hitchin space (2.5c) was also discussed. 
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On the last two curves, the involution acts as z — > —z, with u' invariant. The invariance 
of the curves under the involution means that the equations can be written in terms of v! 
and z 2 . The actual SW curve (quotient by the involution) is thus a cover of the z 2 plane. 
The first curve (2.7a) is invariant under the interchange of the two gauge groups. 

Were we to reverse the transformation, described in the last paragraph, from the curves 
(2.1a)-(2.1c) to the ones in (2.7a)-(2.7c), we would obtain two solutions, due to the fact that 
we squared both sides of an equation in one of the steps above. However, these two solutions 
are related by the involution in cases (b) and (c). In case (a), the two solutions are related 
by e- <-> e-, i.e. P <-> P and so correspond to exchanging the two U(iV) factors. Hence in all 
cases the two solutions are physically equivalent. 

When A' — > 0, the curves (2.7a)-(2.7c) are singular at the roots of P(z), P{z) (or P(—z)), 
and P(z) — P(z) (or P(z) — P(—z)). (The discriminant has double zeros at those points.) 
When A' ^ 0, the surface is deformed such that the first two sets of singular points open up 
into branch cuts, but the singularities at the points z where P(z) = P(z) (or P(z) = P(—z)) 
remain 6 . One (important) exception occurs for the theory with antisymmetric matter, where 
the singularity at z = also opens up into a branch cut in the curve (2.7c). For the theory 
with symmetric matter, z = remains a singular point in the curve (2.7b). 

3 Perturbative approach to the matrix model 

As we discussed in [5] (see also [3]), the M = 2 gauge theory prepotential (in an instanton 
expansion) may be determined using only matrix-model perturbation theory. In this ap- 
proach, one adds to the M = 2 superpotential an additional piece which freezes the moduli 
to a generic, but fixed, point on the Coulomb branch of the M = 2 theory, and breaks the 
M = 2 supersymmetry down to M = 1. After the relevant quantities are computed, the 
extra piece is removed, restoring M = 2 supersymmetry. 

This approach was first explored in [3] where the gauge-coupling matrix was deter- 
mined for U(2), and was extended to U(iV) in [5] where, in particular, a proposal for how to 
determine the relation between the quantum order parameters a$ and their classical counter- 
parts Ci entirely within the context of matrix-model perturbation theory was put forward. 
Using this proposal the prepotential J 7 (a) was calculated to one-instanton level and was 
shown to agree with the well-known result. In [7] the calculation was extended to include 
matter in the fundamental representation, and in [9, 10] to SO/Sp gauge groups. 

In this section, we extend this perturbative matrix-model method to new cases by calcu- 
lating the one-instanton contribution to the M = 2 prepotential in the U(N)x\](N) gauge 
theory with a bifundamental hypermultiplet, and the U(A") gauge theory with one symmetric 
(m) or antisymmetric ([]) hypermultiplet. Besides the additional complication of dealing 
with two-index matter, and the inclusion of diagrams with the topology of RP 2 (for □□ and 
Q), there is one significant modification of the procedure developed in refs. [5, 7]: namely, 
for the models with symmetric or antisymmetric matter, r^- is no longer given simply by the 
second derivative of the free energy [1] (see sec. 3.2.2). Moreover, in the matrix model for 
the U(iV) theory with antisymmetric matter, one is required to expand around a different 

6 This fact is important to get the genus counting to work, cf. sec. 4. 
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vacuum than one would naively have anticipated (this may be related to the results in [22]). 

Previously the one-instanton prepotential for these models has been obtained using M- 
theory methods [31, 32, 33, 34, 35] (see also [37] for another approach). 



3.1 XJ(N)xXJ(N) with a bifundamental hypermultiplet 

Consider the M = 1 U(iV) xU(iV) supersymmetric gauge theory with the following matter 
content: two chiral superfields 0^, 0/ transforming in the adjoint representation of each of 
the two factors of the gauge group, one chiral superfield 6/ transforming in the bifundamen- 
tal representation (□,□), and one chiral superfield bj transforming in the bifundamental 
representation (□, □). The superpotential of this gauge theory is taken to be of the form 7 



W(0, 0, 6, 6) = tr [W(cf)) - W{4>) - 606 + 606], 
where W((f)) and W((f>) are polynomials such that 

N N 

W\z) =aH(z- ej), W\z) = a H(z - ~e 3 ) . 
i=i i=i 



(3.1) 



(3.2) 



The superpotential (3.1) can be viewed as a deformation of an M = 2 theory, which is 
recovered when a — ■> 0. At the end of our calculation we will take this limit thereby obtaining 
results valid in the Af = 2 theory. 

The U(M)xU(M) matrix model associated with this gauge theory [2, 23, 24, 1] has 
partition function 8 



— !— / d$ d$ dB dB expf -— tr[w($) - W($) - B&B + B&B 
vol G J \ 9s 1 



(3.3) 



where $ is an MxM matrix, $ is an MxM matrix, 5 is an MxM matrix, and B is an 
MxM matrix. In the perturbative approach, the matrix integral (3.3) is evaluated about 
the following extremal point of the potential 



/eiV 
e 2 lM 2 








V o 







$0 = 



/MA, o 

e 2 lu 



£n^m n ) 



M 2 








V o 







Bq = B = 0, 



(3.4) 



where = M and E;^i = M. This choice of vacuum breaks the U(M)xU(M) 

symmetry to G = u i M i) x HiLi U(Mi) ■ Writing $ = $ + * and $ = $ + vj>, one finds 
that the off-diagonal fields and (i / j) have vanishing contributions to the quadratic 
part of the action; these fields are zero modes and correspond to gauge degrees of freedom 
[3]. We fix the gauge = = (i ^ j) and introduce Grassmann-odd ghost matrices, 

7 A mass term for the bifundamental fields, mtr(bb), can be introduced by shifting <p — > 4> — m /2 {4> —* 
4> + m/2) and — > — m/2 (e, — > + m/2). 

8 We use capital letters to denote matrix model quantities. 
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exactly as in refs. [3, 5], to which we refer the reader for further details. The bifundamental 
field B (not to be confused with the ghost field in refs. [3, 5]) can be written 



B = 



( B u 

B21 



B\2 

B22 



\B 



Nl 



B 



N2 



Bin \ 
B2N 

Bnn ' 



(3.5) 



where B^ is an Mj x¥; matrix, and similarly for B. Expanding about the vacuum (3.4) 
one finds (in the = ^ = (i ^ j) gauge) 



tr 



B$B - B$B\ = - ejMBjiBij) + ]T tr(5 j ^ ll % - B^^B^) . (3.6) 

Hence, the matrix integral over the quadratic action contains the e^-dependent contribution 
from the bifundamental fields: 



n 



(3.7) 



and the trilinear pieces of (3.6) contribute B^B and B^B vertices to the Feynman diagrams 
(in addition to the vertices considered in ref. [5]). 

We are interested in the planar limit of the matrix model, i.e. the limit in which g s — > 
and Mj, M; — > 00, keeping g s Mi and g s Mi fixed. The connected diagrams of the pertur- 
bative expansion of logZ may be organized, using the standard double-line notation, in a 
topological expansion characterized by the Euler characteristic x °f the surface in which the 
diagram is embedded [38] 

log Z = J297 X F X (S,S) where Si = g s M h Si = g s Mi, (3.8) 

X<2 

where x = 2— 2g with g the genus. In the planar limit, the leading contribution 

F S (S,S) = F x=2 (S,S)=g 2 s log Z 

sphere 



(3.9) 



comes from the connected diagrams that can be drawn on the sphere (x = 2). We will need 
the contributions to F S (S,S) up to cubic order in S and S. The explicit formulae can be 
found in eqs. (B.l) and (B.2) in appendix B. 

To relate the matrix model and its free energy to the Af = 2 U(iV) xU(iV) gauge theory 
broken to IliU(iVj) x IliU(iVj), one introduces, following Dijkgraaf and Vafa, 



W eS (S, S) = - £ N~F s (S, S)-Y: N t ^F s (S, S) 



(3.10) 



where we have dropped terms linear in Si and Si. Since we are examining a generic point 
on the Coulomb branch of the M = 2 theory, which breaks U(iV) xU(iV) to U(1) 27V , we set 
Ni — JVj = 1. Next, one extremizes the effective superpotential with respect to and 5, : 



dW eS (S,S) 



dS,, 



Sj — (Sj)>Sj — (Sj) 



= 0, 



cW eff (S, S) 



dS, 



= 0. 



Sj — {Sj),Sj — {Sj) 



(3.11) 
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The solutions for (Si), (Si) can be evaluated in an expansion in A. The lowest-order contri- 
butions are 



(Si) = ^ N , 



(S t ) = -^A N , 



(3.12) 



where 



N N 

Ti = Y\(e l -~e j ), Ti^Hiii-ej), Ri^Hiei-ej), R t = '[[(e i -e j ) (3.13) 

j=i i =1 j¥=i j¥=i 

and various constants have been absorbed into a redefinition of the cut-off A. In sec. 3.1.2, 
we will also need the next-to-leading-order contributions; these are given in (B.3). 

3.1.1 Relation between a* and e. 

Before computing and the j\f = 2 prepotential, we must determine the relation between 
and the periods a* of the SW differential. In ref. [5], we proposed a definition of a, within the 
context of the perturbation expansion of the matrix model, without referring to the Seiberg- 
Witten curve or differential. As in refs. [5, 7], a* and a, can be determined perturbatively 
via (setting iVj = iVj = 1) 



a, = + 



e» + 



N Q N Q 

E M-9.{K%ds* + E ^ S (tr ^) 52 

;=1 ° b 3 j=l dSj 

N d N d ~ 

E ^9s(tT ^u) S 2 + J^9s(tT ^u) S 2 
3=1 9S 3 j=l db i 



(S),(S) 



(S),(S) 



(3.14) 



where (tr^jj)^ ((tr ^j,)^) is obtained by calculating all connected planar tadpole diagrams 
with an external (^u) leg that can be drawn on a sphere. In addition to the tadpole 
diagrams discussed in ref. [5], there are diagrams with B, B (bifundamental) loops. The 
total contribution to the tadpole quadratic in S and S is 



(tr * ii ) sa 
(tr^u)s2 



ag s 



E 



Si 



Ri^ij Rie 



SiSj 



j Ri hij 



ag s 



q2 qq do. 

_ ^ i _j_ 2 3 — — 



j^i Ri^ij j^i Ri&ij j 



Ri hij 



(3.15) 



ti — e,-, hij = ti — ij and h^ = — e,- = — Inserting these results into 



where e 

eq. (3.14), and evaluating the resulting expression using eq. (3.12), one finds 

0(A 2N ) 



e *+ AA Ml-E r 



1 T T 1 
, _ \p zl il \ " 

d 5 , p2L(,.. 



D2 ^ 

1 j ij 



a. = e "- + A "[|-S^- + iE^-|Ef ]+o(A 2N ). (3-16) 

j^i j ij i j j ij i j ij 
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By using 

' W'(z) 
W'\z) 

one may show that 



E 



Ti 



Ri(z - e;)' 



W'(z) 
W'(z) 



E 



T 

Ri(z - €.i) 



(3.17) 



T 

E 



Ri 



EiT 



-1 



(3.18) 



so that eq. (3.16) may be rewritten 9 

27} 



di = 6i + A 



IV 



1 T- 1 

■^i j^i e ij j ^ij 



Ri 



+ e>(A 2JV ) 



(3.19) 



and similarly for 5j. We cannot yet compare this expression with the SW result obtained in 
[33], because the relation between the roots of W'(z) and the roots e- of P(z) (cf. (2.1a)) 
has not yet been determined. This will be done in sec. 4.1.1. 



3.1.2 Perturbative calculation of t, 



Following Dijkgraaf and Vafa, the gauge coupling matrix is related to the planar free 
energy F s of the matrix model by 



13 2m dSidSj 



d 2 F x 



(S),(S) 



i%3 



2ni dSidS, 



d 2 F K 



(S),(S) 



2ni dSidS, 



(S),(S) 



(3.20) 



We may calculate these expressions perturbatively using eqs. (B.l), (B.2), and (B.3), and 
finally, use eq. (3.16) to re-express the entire expression in terms of rather than e^. The 
resulting perturbative and one-instanton contributions to the gauge coupling matrix are 
given in eqs. (B.5)-(B.7) in appendix B. One may verify that (B.5)-(B.7) can be written as 



d 2 J r (a ) a) 
ddjddi 



d 2 T{a, a) 
ddidcij 



d^Fja, a) 
ddiddj 



(3.21) 



with (up to a quadratic polynomial) 



2TriJ r (a,a) = -i ^ ^(a* - a^) 2 log 



+iE( a »-%) 2l °g 



a,- — a. 



A 



iEE(^- fi i) 2l °g 

i j^i 



OLi — dj 



A 



a,; - a,, 



A 



+ A"E 



R) 



+ 0(A 2N ) . (3.22) 



which agrees perfectly with (version 2 of) ref. [33]. 

9 In previous work [5, 7] we used similar identities at intermediate stages of the calculations. However, in 
the calculation of the J\f — 2 prepotential, it is more efficient to work with the expressions that come naturally 
out of the matrix-model calculation. To compare with results obtained using M-theory at intermediate stages, 
identities generally have to be used (as we did to obtain (3.19)). 
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3.2 U(N) with m or g 

Consider the Af = 1 U(iV) supersymmetric gauge theory with one chiral superfield trans- 
forming in the adjoint representation of the gauge group, one chiral superfield Xij trans- 
forming in either the symmetric (m) or the antisymmetric (Q) representation, and one 
chiral superfield transforming in the conjugate representation. We treat the cases of the 
symmetric and antisymmetric representations simultaneously by assuming that x, x satisfy 
x T = f3x and x T = (3x, where /3 — 1 for the symmetric representation and f3 — — 1 for the 
antisymmetric representation. The superpotential of the gauge theory is taken to be of the 
form 10 

W((j),x,x) =tr[W(0) -x<j>x\ , (3.23) 

where W((f>) is a polynomial such that W'(z) = otY[f = i(z — ef). This superpotential can be 
viewed as a deformation of an Af = 2 theory, which is recovered when a — > 0, restoring 
Af = 2 supersymmetry. 

As discussed in [25] there are several classical ground states of the superpotential (3.23). 
One such ground state is = e^l and x = x = 0. Another one is = 0, x = E and 
x = W'(ti)E~ l , where E — t for m and E = J for [], where J is the usual Sp-unit. 
There are also additional ground states as discussed in [25] , but these will play no role in our 
discussion. (Similar extra vacua are also present [39, 23] in the U(iV) xU(iV) theory discussed 
above.) A more general vacuum is obtained by combining the above possibilities. In a block- 
diagonal basis, one ground state is = diag(0^ o , eil^, . . . , efcljv fc ), with N = N + Yh=i Ni 
(where N is even for Q) and x and x have vanishing entries except for the N xN blocks 
xqo = E, x 00 = W f {0)E~\ Such a vacuum breaks XJ(N) down to [25] SO(A" ) x Hi U(A r i) for 
□ or Sp(iVo)xn j U(iV l ) for g. 

We want to freeze the M = 2 moduli to a generic, but fixed, point on the Coulomb branch 
of the Af = 2 theory. This is accomplished by breaking U(iV) down to XJ(1) N , i.e. choosing 
Ni = 1 and N = 0. 

The U(M) matrix model associated with this gauge theory [25, 1] has partition function 11 

1 f .„ ~ (I 



J d$dXdX ex pf-— tr[w($) -X<Z>X]\ , (3.24) 



vol G J \ g 

where X T = f3X and X T = (3X. In the perturbative approach, the matrix integral (3.24) is 
evaluated about an extremal point of W(Q,X, X). 

Based on previous experience, it would seem natural to expand around a matrix-model 
vacuum similar to the gauge theory vacuum but with Ni — 1 and No — replaced by Mj 
and M such that Mj ^ and M = 0. This will indeed turn out to be the right procedure 
for the U(A r ) + □□ theory. However, as we will see, it is not the right procedure for the 
U(iV) + g theory. Instead, for this theory, we will take Mo 7^ (in fact we will take the limit 
M — > 00, with g s M finite), even though N = 0. We do not have an a priori reason for 



10 



A mass term mtr(xx) for the matter hypcrmultiplct can be introduced by shifting <p — > <p — m and 

e« — * £i — m. 

11 As in the previous section, we use capital letters to denote matrix model quantities. All matrix indices 
run over M values. 
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making this choice of vacuum. If one does not include the extra Mo x Mq block for the matrix 
model corresponding to the U(iV) + Q theory, one still gets an (apparently) self-consistent 
result, but one which does not agree with the prepotential, SW curve, or SW differential 
derived from M-theory [29, 31, 37] (which have passed several consistency tests). Only if one 
includes the extra block does one get a result that is in agreement with previous results. 12 
We will decompose all matrices T as 



T 



/ T~oo T i 
Tio T n 



Yin 



(3.25) 



V T jvn T jvi ' • • Y Any / 



where is an Mj xMj matrix, Tj is an MjXM matrix, T o is an M xM matrix (where 
M is even for Q). Throughout we use i, j = 1, . . . , N, displaying the index-0 terms explicitly. 
We evaluate the matrix integral (3.24) about the following extremal point of W(<&, X, X) 



and 



(3.26) 







/ Mo 




^ 








$ = 


eilMi 







J 








v 


••• £N^M n J 








■•• 


o\ 






( E ~ x 


• 


■ °^ 





••• 





i X Q — 


W'(0) 





• 


• 


U 


••• 


o) 






I 


• 


• o) 



(3.27) 



where M + YliLi M = M, and as before E = l Mo for CD, and E = J for where J is 
the antisymmetric Sp(M ) unit. This choice of vacuum breaks the U(M) symmetry to [25] 
G = SO(M )x nili U(M-) for m, and G = Sp(M )x Uti U(Mi) for Q. 
The potential in (3.24) is invariant under the gauge symmetry 



<5<t> = [£,<£>], 5X = £X + xe, 5X = -eX-X£. 



(3.28) 



Writing $ = $ + * and X = X + Y, we fix the gauge ^ = (ij^j), = 0, = 0, 
F o = 0. Following [25] we use the BRST approach and introduce the gauge- fixing fermion 

6 = Vtr(C ,<!>,/) + "> " ! \:(C ', ( i<I> n , ) + " 1 i-(C' n , <[>,■(,) -f 1 r(( ' nll X m — E}) (3.29) 



5>(C^) + E tr (a <fo,) + £ tr (<7<A) +tr(C 00 [X ( 



where the C's are Grassmann-odd. The relevant 

s$ = [C,$], sX = CX 

12 Added in v4: Although the choice of vacuum appears 
it is quite likely that the vacuum is uniquely determined 
Abel's theorem (method II in section 4.2.2 of this paper), 
periods of T(z) are integer- valued [40]. See refs. [41] and 
the vacua in related theories. (The M x M block must 
in the algebraic curve around = 0.) 



BRST transformations are [25] 
+ XC T , sC = D 



(3.30) 



ad hoc within the purely perturbative framework, 
through the exact determination of the curve via 
or equivalently, through the condition that all the 
[42] for the use of the latter criterion to determine 
be included whenever an additional cut opens up 
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where C and D are Grassmann-odd matrices. The gauge-fixing action S g { — s G becomes, 
after using (3.30) and integrating out the .D's which act as Lagrange multipliers implementing 
the gauge choice 13 

ijtj i i 

+ tr(C 00 [C 00 E + EC^\) + 5>(q,o[a,^o + *o*Q£]) (3.31) 



Expanding about the vacuum (3.26), (3.27) using the above gauge one finds (in addition 
to the quadratic terms in eq. (3.31)) the quadratic part of the action 



N 



tr[W($) - X$X] 



\a J2 RM^l) + |«i?otr(^ ) - tr(F o^oo^) (3.32) 

i=l 

N N 

- eMY 0i Y i0 ) - e^tr (y^) - ^(e, + e,)tr(Y^) + • • • 



i=i 



where 



Ri 



IK* 



i<j 



(3.33) 



From this we see [25] that the antisymmetric matrix Y 00 acts as a Lagrange multiplier 
implementing the constraint ^ oE + /3(^ 00 E) T = ^ 00 E + E^ = 0, i.e. ^ o e so(M ) for 
m and ^oo e sp(M ) for Q. 

The matrix integral over the quadratic action contains the e^-dependent contributions 



Rq 



n 



i \ i M » 2 / 1 \ \ M i( M i+p) 



R; 



ei) 



M Mi 



n 



MiMj 



C{ ~\~ Cj 



(3.34) 

The following terms in the action (in addition to those considered in ref. [5] and the trilinear 
terms in eq. (3.31)) contribute cubic vertices to the Feynman diagrams: 

N 

E [ tr (£o*oo*oi) + tr (Yo^Fo) + tr (%y u Y u )] - ^trfe**^ + Y^ 33 Y 3l ) . (3.35) 



i=l 



We are interested in the planar limit of the matrix model, i.e. the limit in which g s — > and 
Mj, M — > oo, keeping g s M^ g s M fixed. The connected diagrams of the perturbative ex- 
pansion of log Z may be organized, using the standard double-line notation, in a topological 
expansion characterized by the Euler characteristic x °f the surface in which the diagram is 
embedded [38] 

log Z=J2g;*F x (S,S) (3.36) 

where x — 2— 2g— q with (7 the genus (number of handles) and q the number of crosscaps. In 
the Feynman diagrams, we generally replace g s Mi by Si and g s M by So, but for the inner 



13 



In refs. [3, 5] the notation B, C was used for the ghost fields instead of C and C. 
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index-loop of an X^, X^ loop we write g s Mi = Si since the arrow on the inner index-loop 
runs parallel to the outer index-loop, opposite to the direction in which it would run for 
the adjoint representation This will be important when (but not until) we calculate in 
sec. 3.2.2. 

In the planar limit, the leading contribution to the matrix integral comes from the planar 
diagrams that can be drawn on the sphere (% = 2), 



F s (S,S) = F x=2 (S,S)=g 2 s \ogZ 



(3.37) 

sphere 



The sub leading contribution comes from planar diagrams that can be drawn on PJP 2 , which 
is a sphere with one cross-cap inserted (x — 1) 



F rp (S) = F x=1 (S)=g s \ogZ 

EP 2 

To evaluate the (1/vol G) prefactor in (3.24) we need (see e.g. [43, 44, 25]) 



(3.38) 



logvol(SO(M )) = ^ M 2 log M + ±M log M + ... 
logvol(Sp(M )) = 4 M 2 log M -±M log M + ... (3.39) 
logvol(U(M i )) = -|M 2 logM i + ... 

These results together with the integration over the quadratic fields of the matrix-model 
partition function yields (up to an e^-independent quadratic monomial in the S"s) 

N N / <7- \ N /p-\ 



\ E log g ) - \ e E sft log (f ) + E SoSi io g (| 



n>3 



where e^- = — and = e« + e^. The term F^ n \S, S) is an nth order polynomial in Si 
and Si arising from planar loop diagrams built from the interaction vertices. We will also 
need the cubic contribution F^ 3 \S,S); the explicit expression can be found in eq. (B.8) in 
appendix B. 

Next we turn to the sub-leading contributions to the free energy. From eqs. (3.31), (3.32), 
(3.39) one finds (up to an ej-mdependent part linear in Si) 14 



Fr P (S) = -^E^aog(|)-^olog(-^) + E^ n) (^) (3.41) 



We note that the sphere and MP 2 contributions to the free energy obey the relation [45, 44, 46] 

^(S) = -l^-F s{ S,S). 
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where F^(S) is an nth order polynomial in S, arising from planar diagrams, built from the 
interaction vertices, that can be drawn on RP 2 . We will need the quadratic contribution 
F r (2) (S); the explicit expression is given in (B.9). 

To relate the matrix model to the M = 2 U(iV) gauge theory broken to IliU(iVj), one 
sets S — S in the matrix-model free energy [1] and introduces [45, 44, 46, 25, 1, 19] 

We S (S) = -^Ni—F s (S,S)-N —F s (S,S)-4F rp (S) (3.42) 

where we have dropped terms linear in Si. Since we are examining a generic point on the 
Coulomb branch of the M = 2 theory, which breaks U(iV) to XJ(1) N , we set iVj = 1 and 
Nq = 0. Next, one extremizes the effective superpotential to obtain (Si) and (So): 



dW eS (S) 



dSi 



0j dW eS (S) 



S J -=<S J ->,So=<So> ' «9S, 



= 0. (3.43) 

Sj=(Sj),So=(So) 



The solution for (Si) and (S ) can be evaluated in an expansion in A. The lowest-order 
contribution is 

( Si ) = ^A N -W, G t = ef YKei + ej), ^ = 11^-^) (3-44) 



and 



(S'o) = -a(-l) N A^ N+2 W"(0) (^Mj (3.45) 



and constants have been absorbed into a redefinition of the cut-off A. The next-to-leading 
contribution, which will be needed in sec. 3.2.2, is given in (B.10). 

Now for (S'o) and (Si) are both 0(A N+2 ) and therefore both need to be included in 
a perturbative computation to this order. For l l l , however, (So) is inversely proportional 
to A, which seems to indicate some inconsistency in perturbing about the vacuum (3.26), 
(3.27) when (5 — 1. Therefore, in the □□ case, we will simply expand around the vacuum 
with Mq = instead; equivalently, we will use 

(S )=2a5^ 1 A N + 2 (-l) N Y,- (3-46) 



in all expressions below. 



3.2.1 Relation between a, and e, 



Before computing r^- and the M = 2 prepotential, we must determine the relation between 
Oj and ti. As in refs. [5, 7], a, can be determined perturbatively via (setting iVj = 1 and 
A r o = 0) 

N g 

di = ei+ J]— ff,(trf ii ) S 2 + 4(tr* ji ) Bpa (3.47) 
J =1 j J {S) 
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where {tr^/u) S 2 is obtained by calculating all connected planar tadpole diagrams with an 
external leg that can be drawn on a sphere, and (tr \I/ 33)^2 is obtained by computing 
all connected planar tadpole diagrams with an external leg that can be drawn on HP 2 . 
(The factor of 4 in eq. (3.47) arises from the corresponding factor in W e g (3.42), using the 
arguments in ref. [5].) The total contribution to the tadpole on the sphere quadratic in Si is 



(tr^) 5 2 



1 



ag s 



E 



j Ri9ij 



RiC{ 



(3.48) 



and contribution to the tadpole on HP 2 linear in Si is 

P Si 



(tr *ii)np2 



(3.49) 



Inserting these results into eq. (3.47), and evaluating the resulting expression using eqs. (3.44) 
and (3.46), one finds 



2 x Gj 

E 



_ J_ v Gj 



Ri j-£i -^j^ij -^i j -^jdij ^'i -^"> 



Gi 1 2[3Gi 
j ■ l 'J 



3 e i 



(3.50) 



To compare these results with those found using SW theory, we must consider the symmetric 
and antisymmetric cases separately. 

(i) u(iv) + m (P = 1) 

Consider the function 



N 



- z 2 - 2a 1 z - 2al 



(3.51) 



where 0\ = X^Li e ii an d the last three terms remove the no n- negative powers of the Laurent 
expansion of g(z). The function g(z) has only simple poles at z = ti and so can be written 
in terms of its residues as, cf. (3.44), 



= E 



Gi 



Ri(z — e 

Using eqs. (3.51) and (3.52), one may show that 

1 



E Gj 

j^i Rj e ij 

E Gj 

j Rj9ij 



Gi 
Ri 



-E 



1 2 



= e\-2o x ei + 2o\, 



e 2 - 2a 1 e i - 2a\ , 



(3.52) 



(3.53) 
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so that using eq. (3.50) we get 

A N-2 



di = e; + 



Ri 



— 3e 2 — 2a"! ej — < oo\ 



it,- 



2 1 2 



+0{A 



2iV-4\ 



(3.54) 



After determining, in sec. 4.2.1, the relation between the roots of W'(z) and the roots 
of P(z) in (2.1b), we will be able to compare this result with that obtained from SW theory. 

(ii) U(AT) +0(0= -1) 
Consider the function 

■W'(-z) 



h(z) = (-1) 



N 



W'(z) 



Now h(z)/z 2 has a double pole at z — 0, which we may remove by writing 



- h(0) - zh'(0) 



\ A' 



N 



W'(z) 



-l-2zJ2 



i=i e * 



The function if(-z) has only simple poles at z = e« and so may be written 

H (z) = V — ^ - . 

Using eqs. (3.56) and (3.57), one may show that 



Ri 



\N 



so that using eq. (3.50) we get 
3(-) N A N+2 



1 1 2 

E-+E--- 

j 9ij &i 



l + 2e^- 



1 



2e,Ef 



e,; - 



^ 2 ^ , 1 2 

7 - fl'ij 6j 



+ 0(A 2JV+4 



(3.55) 



(3.56) 



(3.57) 



(3.58) 



(3.59) 



This equation, obtained entirely using matrix model methods, precisely agrees (after letting 
A = —A', cf. eq. (4.43)), with eq. (4.2) in ref. [31], obtained using the Seiberg-Witten 
procedure, because, as we will see in sec. 4.2.1, the roots e, of W'(z) and the roots of P(z) 
in (2.1c) coincide to this order in A. 



3.2.2 Perturbative calculation of 

For these models, in contrast to models containing only adjoint and fundamental represen- 
tations [15, 13] or bifundamental representations (sec. 3.1), the gauge coupling matrix r^- is 
not given by the second derivative of the planar free energy F s of the matrix model: 



d 2 F !i 



Sk — {Sk) 



(3.60) 
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Nevertheless, using a diagrammatic argument, a perturbative prescription for can be 
given [1] 

1 ( d d \ ( d d 



13 



F S (S,S) 



Sk—Sk — {Sk) 



(3.61) 



J 3 ^ /LJ 3 / 

(It is unclear what, if any, the physical meaning of derivatives w.r.t. So is-) 

In appendix B we evaluate the expression (3.61) perturbatively up to one-instanton order, 
expressing the result in terms of a, rather than e^. The result can be written as = 
d 2 J r (a)/da i daj with (up to a quadratic polynomial) 



i j^i 



log 



A 



+ IEE(^ + %) 2l °g 

i j^i 

25/3-1 



Hi ai 



cii + a. 



+ 0(A 2N ~ 4/3 ) . 



(3.62) 



Later (cf. eqs. (4.37, 4.43)) we will see that A = —A', where A' is the quantum scale used in 
the M-theory curves (2. lb, 2.1c). Taking this into account, eq. (3.62) precisely agrees with 
the calculations of the prepotential in refs. [31, 32] which utilize the SW curves for these 
theories derived from M-theory [29]. 



4 Cubic matrix- mo del curves 

In this section we study the algebraic curves that arise from the planar solution of the 
matrix models and show how to obtain from these the SW curves of the M = 2 gauge 
theories discussed in sec. 2. 



4.1 XJ(M)xU(M) matrix model with bifundamental matter 

The large-M planar solution of the U(M)xU(M) matrix model described in sec. 3.1 was 
discussed in [24, 1] (several of the results can also be found in [47, 2, 23]). In this approach, 
one defines the resolvents 



(4.1) 



n=0 x \^ ^V' n=0 

where matrix-model expectation values are defined via 

(o($, $, B, B)j = i J d$ d$ dB dB $, B, B) e ~j;H w ^-^~^ B + B ^} . (4.2) 
It may be shown that 

u 1 (z) = -u(z)+W'(z), u 2 {z) = u(z) - Cj(z) , u 3 {z)=u>(z)+W'{z) (4.3) 

(where oj{z) is the leading (sphere) part of the resolvent) are the values of a variable u on 
the three sheets of a Riemann surface defined by 15 

(u - W'{z)) u{u- W\z)) = n(z) u - h(z) . (4.4) 

15 This equation may be obtained e.g. from eq. (2.22) in ref. [1] by redefining u — > —u + 3 (W'(z) + W'{z)) 
and setting h{z) = si(z) + \{W'{z) + W'(z))n(z) 
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The coefficients of the cubic curve (4.4) are given by [24, 1] 
ri(z) = -g s (tr(^ 

and 

tiW = -&wc*)w V;_ $ - l)+g«^)(tr( 1 - j) 

+ M w ' iz i - r* V w ) ) - 4 r p /( t : r*^ ) ) (4 - 6) 

from which one can see, using eq. (3.2), that r\(z) and are polynomials of degree at most 
N—l and 2N—1, respectively, whose coefficients depend on the vevs (tr($ fe )) and (tr($ fc )) 
with k < 2N—1. At this point these vevs, and therefore r\(z) and t 1 (z), are undetermined. 
We would now like to connect the above general curve to the cubic Seiberg-Witten curve 
(2.7a) for the Af = 2 theory. We will discuss two methods. 



\\ / (W'(z)-W'(&)\\ 



_d_ ( W'{z)-W'{$) 
d$ I z - $ 



^ 2 (tr 



_d_ - 
d$ V z - $ 



4.1.1 Method I: perturbative determination of the curve 

The planar solution to the matrix model yields a curve (4.4) dependent on arbitrary poly- 
nomials ri(z) and t\{z). As we saw above, the coefficients of these polynomials depend on 
(tr($ fc )) and (tr($ fc )), which at this stage are arbitrary. An additional condition is necessary 
to fix these polynomials, namely, the extremization of W e s- This will determine (tr($ fc )) and 
(tr(<3> fc )), and thus lead to a specific form of the cubic curve. Only then can the matrix-model 
curve be compared with the Seiberg-Witten curve obtained from M-theory (see section 2). 

One method of using the extremization of W e s to determine the curve employs Abel's 
theorem, and was described in section 7 of ref. [7] (see also [4, 9]). This approach will be 
discussed below in subsection 4.1.2. 

However, the method using Abel's theorem is difficult to apply in some cases, so in this 
section we will present an alternative approach that is more straightforward to implement. 
This method is to evaluate (tr($™)) perturbatively in powers of A, and use the result to 
determine r±(z) and ti(z), and therefore the form of the cubic curve, order-by- order in 
perturbation theory. Although this method does not yield the exact form of the curve, it is 
a quick and efficient way of determining the form of the curve to lowest order in A. 

Expanding $ = $ + where $ is given by (3.4) one easily sees that to lowest order in 
perturbation theory the matrix model expectation values g s (tr($ n )) are given by J2i{Si) e i- 
Thus, writing 

w'iz) - w'm ^ , , _ n , . 



* - $ n=0 



and using W'{ei) = 0, we have 



i <f> / / ^— ' \ /j" \ \ ii / / \ <■ i l 4 11 \ / % \ j / ^ 
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where the second equality only holds to lowest order. If one is only interested in the lowest- 
order contribution, one can drop the last four terms in ti(z) (4.6). (The terms on the second 
line of (4.6) are double-trace terms; since they contain products of at least two S"s, they are 
at least second order in A N . The terms on the last line of (4.6) vanish to lowest order since 
W'{ei) = W'{ei) = 0.) Using (4.8), one finds (to lowest order) 



ri (z) = -r(,)E^ + rwE {Si) 



h(z) 



-W'(z)W'(z) ]T + W\z)W'{z) ]T -^z 



(4.9) 



Using eq. (3.12) we find, again to lowest order, 



Y Ri{z - et) ' 

Inserting eq. (4.10) into eq. (4.9) and using eq. (3.17), one obtains 



ri ( z ) = + O(A 



2N\ 



ti(z) = -aA N W'(z) - W'(z) + 0{A 2N ) 



hence the matrix model curve is (to first order in A N ) 

(u - W\z)) u(u- W\z)) = aA N [W(z) - W\z)] 2 + 0(A 2N ) 
This curve is identical to the (transformed) M-theory curve (2.7a) provided that 

W'{z) = a (P(z) - 3A' N ) , 



u = au 



W\z) = a(P(z)-3A' N ), A = A' 



(4.10) 



(4.11) 



(4.12) 



(4.13) 



In summary, the matrix model curve (4.4), together with the extremization of W e g, which 
gives (3.12) and therefore (4.11), leads to the U(iV)xU(iV) SW curve (2.7a). 

The relation (4.13) implies that the roots of the polynomial P(z) = YiiLi( z — e-) in the 
SW curve (2.1a) and the roots of the derivative of the matrix model potential W'(z) = 
a HiLi( z — e i) are equivalent classically, but differ by 



O A N 



(4.14) 



to first order in A^, and analogously for t { and e\. This just amounts to a redefinition of 
the moduli and h L . In sec. 3.1 we determined the relation (3.19) between the SW periods 
<2j and ej. Combining (3.19) and (4.14) allows us to write 



* = < + A /JV f 



1 1 



+ 0(A' 2 ") 



(4.15) 
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(and a similar relation between on and eQ. This results precisely agrees with eq. (17) of 
ref. [33], obtained using the Seiberg-Witten procedure. 

The fact that the first-order curve (4.12) precisely agrees with the M-theory curve (2.7a) 
(which is believed to be the exact answer) points to the existence of a non-renormalization 
theorem in the matrix model (which we have not proven). As shown in appendix A, a similar 
result holds for U(iV) with Nf fundamentals when Nf < N (but not when Nf > N). Thus 
in some cases (but not always) the perturbative method described above actually gives exact 
results. 

4.1.2 Method II: exact determination of the curve via Abel's theorem 

In this section, we will follow the strategy of [7] and discuss, using the saddle-point solution, 
the condition on the matrix-model curve imposed by extremizing W e s. 

The cubic curve (4.4) with the right-hand side set to zero is a singular curve with sin- 
gularities at z = ti (the roots of W'(z)), at z = (the roots of W'(z)), and at the roots of 
W'(z) — W'(z). Turning on the right-hand side generically deforms the curve into a three- 
sheeted Riemann surface with (square-root) branch cuts between sheets one and two located 
near e^, branch cuts between sheets two and three located near e^, and branch cuts between 
sheets one and three located near the roots of W'(z) — W'(z) (see, e.g., [23] for a picture 
of the cut structure of this curve). This generic Riemann surface has genus 3 A — 2. If, 
however, the last described set of cuts does not open up, the curve remains singular, having 
(geometric) genus 2N — 2. (This is in fact the case for the SW curve (2.7a), agreeing with 
the fact that the M = 2 moduli space is 2(N — l)-dimensional.) 

To impose the extremization of W c s on the matrix model curve (4.4), we begin by ex- 
pressing the leading (sphere) part of the free-energy of the matrix model in an eigenvalue 
basis as (cf. [23]) 

F s = JdX dX' [p(A) P (Y) log(A - A') + p(A)p(Y) log(A - A') - p(X)p(X') log(A - A')] 

- JdX[p(X)W(X)-p(X)W(X)} (4.16) 

where p(A) and p(A) are the densities of eigenvalues 

p(A) =g s J2 S(X - A,) , p(X) =g s J2 S(X - K) (4.17) 

i i 

(normalized as / dA p(X) = g s M = S and / dXp(X) = g s M = S), which are related to the 
resolvents (4.1) via 

uj(z) = [ dA^- , Q(z) = f dX^- . (4.18) 

J z — X J z — X 

Next we define 16 

Si = -^-.i u{z)dz, S i = ^-£u(z)dz (4.19) 

2m JAi Zni JAi 

16 In general there are also corresponding S's for the cuts connecting sheets one and three (see [23] for a 
discussion). However, these will not affect our discussion so we will suppress them. 
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where Ai and Ai denote contours around the branch cuts near e« and on sheets one and 
three respectively. Using (4.3) and (4.18) one can show that Si and Si are the integrated 
densities of eigenvalues along the cuts near and e«. (Thus the definition (4.19) is consistent 
with the perturbative definition Si = g s Mi, Si = g s Mi.) Variations in Si and Si can be 
implemented [48, 7] by varying the densities Sp(X) = 5 Si S(X — e«) and Sp(X) = 5 Si 5(X — e^), 
with Ci and denoting any point along the branch cuts (cf. [40] for an alternative approach). 
Specifically, (up to terms which will not affect our discussion; see [40] for a discussion of such 
terms) 



m 

dS t 



-W(e t ) + J dX[2p(X) log(A - a) - p(X) log(A - e,)] 
= f dzW\z) + [dX -2p{X) f dz-^— + p{X) f dz— 
= f dz[W'{z) -2oo{z)+u{z)}= f dz[ Ul {z) -u 2 {z)}= j h u (4.20) 

J d Jet Jh 

where we have used (4.3), and the last expression is interpreted as the integral of u from 
I 2 , infinity on the second sheet, to I±, infinity on the first sheet, along a contour that passes 
through the cut near e^. Similarly, 

dF r 1 ^ 

—-lr- — J u (4.21) 

dSi Jh 

where the integral is taken along a contour that passes through the cut near e^. The results 
(4.20), (4.21) were also obtained in [23]. 

Next, we wish to extremize the effective superpotential (3.10) 

N dF N dF 

(setting Ni = Ni = 1) for the U(iV) xU(A r ) theory broken down to U(1) 2JV . This is accom- 
plished by taking derivatives of W e s w.r.t. the Si's and Sj's. However, in analogy with [4, 7], 
one may change variables and instead vary w.r.t. the coefficients of the arbitrary polyno- 
mials r\(z) and t\(z) in the matrix model curve (4.4). From (4.4), one may check that the 
derivatives of u w.r.t. the coefficients of ri(z) and t\(z) are holomorphic differentials on the 
Riemann surface 17 . We can change basis to the canonical basis of holomorphic differentials, 
(k, dual to the homology basis, so that the conditions arising from extremizing W e s may be 
written as (see ref. [7] for further details) 



rh 

N Cfc = (modulo period lattice) . (4.23) 



This condition implies, by Abel's theorem, the existence of a function with an iVth order 
pole at Ii, an iVth order zero at I3, and regular everywhere else. For a generic choice of 
ri(z) and ti(z) such a function will not exist (by the Weierstrass gap theorem). Thus only 
in very special circumstances can such a function exist. 



17 Actually not all these differentials arc holomorphic; we only consider the holomorphic ones. 
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Let us first show that the problem has a solution. Consider eq. (4.4), with r\(z) = 
and 18 ti(z) oc (W'(z) — W'(z) 2 , i.e. a curve of the form: 

[u - W'{z)) u{u- W'{z)) = const x (W'(z) - W\z) f . (4.24) 

One can infer the asymptotic behavior of u at 7j, infinity on each of the three sheets 

h ■ u = W\z) + Oiz' 1 ) , I 2 : u = 0(z~ x ) , h ■ u = W\z) + 0{z~ l ) (4.25) 

(Alternatively, this can be deduced from (4.3) together with the asymptotic behavior of the 
resolvents.) 

Now consider the function, defined on the curve (4.24), 

This function has the right asymptotic properties to satisfy Abel's theorem, but potentially 
has poles and zeros at finite values z . The denominator vanishes at any point z at which 
u(zo) = W'(zo) on one of the sheets. However, using (4.24), one sees that the function is 
actually regular at these (singular) points. Similarly the potential zeros at the points where 
u(zo) = W'(zo) are absent. Thus the function (4.26) defined on the curve (4.24) has the 
divisor implied by (4.23) via Abel's theorem. 

Assuming that the solution is unique, we find that the matrix model implies a Seiberg- 
Witten curve of the form (4.24) (with no terms higher order in A). This precisely agrees with 
the M-theory result (2.7a) after the redefinitions (4.13). Note that after these redefinitions, 
the function (4.26) is proportional to the variable y (2.6a) appearing in the M-theory curve 
(2.1a). We have not been able to show uniqueness. 



4.2 U(M) matrix model with symmetric or antisymmetric matter 

The large-M planar solution of the U(M) matrix models described in sec. 3.2, was discussed 
in refs. [25, 1] (these models are a slight modification of the 0(1) model described in refs. [49]). 
In this approach, one defines the resolvent 

u(z) = g.(tr (j^)) = 9s £ ^ 1 (tr($™)) (4.27) 

where matrix-model expectation values are defined via 

(c?($, x, x)) = | J d$ dx dx e>($, x, x) e -^r[w(<E>)-x$x] _ (4 28) 

Then it may be shown that 

u 1 (z) = -uj (z) + W > (z), u 2 (z) = uj (z) + ojo(-z) , u 3 (z) = -u (-z) + W'(-z), (4.29) 

18 If one starts with an arbitrary t\{z) instead, the requirement that the function in eq. (4.26) below should 
have the desired properties implies that t\(z) has to be of this form. 
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where ujq(z) is the leading (sphere) part of the resolvent, are the values of a variable u on 
the three sheets of a Riemann surface defined by 19 

(u - W'{z)) u(u- W'(-z)) = n(z) u - h(z) (4.30) 

The coefficients of the cubic curve (4.30) are given by [25, 1] 

= -*(*( ))- s 'H - J) (4 - 3i) 

and 

" d (W'(z) - W'(<Z>Y ' 



. W'iz) - 11 '(<!>) , 



+gs( tr( ^:f^ ^ 



d$ \ z — <& 
+ {z^ -z) (4.32) 



from which one sees that r\(z) and ti(z) are even polynomials of degree at most N—l and 
2N—2, respectively, whose coefficients depend on the vevs (tr($ fe )) with k < 2N—1. 

4.2.1 Method I: perturbative determination of the curve 

As in section (4.1.1), we may evaluate the polynomials r\ and t\ perturbatively in A, and 
use the result to determine the curve order-by-order in perturbation theory. Expanding 
$ = $o + where <3> is given by (3.26), one sees that, to lowest order in perturbation 
theory, the matrix model expectation value g s (tr($™)) is given by J2i(Si) e" + (So)<5n,o- Hence, 
similar to eq. (4.8), we find 

ri ( z ) = -W'(z)J2^--W'(z)^- + (z^-z). 
i z — ej z 

tl (z) = -W'(-z)W'(z)J2^- -W'(-z)W'{z)^- + (z^ -z). (4.33) 

j z — ti z 

Using the lowest-order perturbative results (3.44), (3.46) we have 

£ J&L = aA^ £ ^T^T , <*o> = 2a 5^A^(-lf £ \ . (4.34) 

i 1 i ^ V 1 J j 3 

To evaluate J2i[Gi/Ri(z — e;)], we must consider the U(iV) + □□ and U(iV) + cases 
separately. 

(i) \j(n) + m 09 = i) 

Using eqs. (3.51) and (3.52), we find (to lowest order) 

ri (2) = -aA^W"^)^*) + (*->-*) (4.35) 
= aA^ 2 (z 2 - (-lfz 2 + 2a 2 ) [W\z) + W\-z)\ + 2aA N ~ 2 a 1 z [W'(z) - W'(-z)} 



19 This may be obtained, e.g., from eq. (3.20) in ref. [1], by redefining u — > — u + \{W'{z) + W'{—z)) and 



setting h{z) = Sl (z) + \{W'{z) + W'(-z))n(z). 
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and (also to lowest order) 

h(z) = -aA N - 2 W'(-z)W'(z)g(z) + (z^ -z) (4.36) 
= -aA^-^-l)^^ 2 [W\z)-W\-z)] 2 + 2aA N ~ 2 (z 2 -(-l) N z 2 +2af) W\z)W\-z) 

The bottom lines of these equations make manifest that r\(z) and t\(z) are (even) polyno- 
mials, and the top lines show that they are of degree N — 1 and 2N — 2 respectively. 

The cubic equation (4.30), with r±(z) and ti(z) as above, may be considerably simplified 
by defining 



W'{z) 



u 
A 



a 



a 



P(z)+A 



N-2 



-3z 2 — 2o\Z — 6a 2 



0(A 



27V+4\ 



(4.37) 



u' - 2A N - 2 {z 2 - {-l) N z 2 + 2<r 1 2 )l + 0(A 2N+4 ) , 



-A'. 



Then 

(v! - P{z) + 3A ,7V " V) v! («' - P(-z) + 3A ,N ~ 2 z 2 ) = A' N - 2 z 2 {P{z) - P{-z)f + 0(A' 2N ~ 4 ) 

(4.38) 

in agreement with the transformed M-theory curve (2.7b). 

The relation (4.37) implies that the roots of the polynomial P(z) = Y\iLi(x— e'j) in the SW 
curve (2.1b) and those of the derivative of the matrix-model potential W'(z) = Ylf = i(x — e^) 
are equivalent classically, and are related by 



e,- = e { + 



A 



N-2 



Ri 



-3e 2 - 2a iei - 6a 2 ) + 0{A 2N ^) 



(4.39) 



at the one-instanton level. Combining this result with the relation (3.54) between a« and e^, 
we find 



e> + (-) V~ 2 § 

it; 



E2 ^ — > 1 2 



+ 0(A 



I 2N-A\ 



(4.40) 



jfr ^3 j 9ij 

This equation precisely agrees with eq. (26) of ref. [32], obtained using the Seiberg-Witten 
procedure. 

(ii) u(Ao + go3 = -i) 

Using eqs. (3.56) and (3.57), one finds (to lowest order) 



n(z) 



-aA N+2 W\z) 



H(z) + 2(-lfWy 

Z „• Cj 



+ 



= 



(4.41) 



and 



t^z) = -aA N+2 W'{-z)W'(z) 



H(z) + 2(-ifWy 

Z ■ K,; 



+ (z - -z) 



W'(z) - W'(-z) 



1 2 



(4.42) 
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Defining 

W\z) = aP(z) , u = au', A = -A' , (4.43) 

equation (4.30) may be rewritten as 

A/iV+2 

(u 1 - P(z)) u (v! - P(-z)) = —^—(P(z) - P(-z)) 2 + OiM 2N+A ) (4.44) 

z 1 

which is in agreement with the (transformed) M-theory curve (2.7c). 

Note that (4.43) implies that the roots of the polynomial P(z) = ]\f = i(x — e-) in the SW 
curve and those of the derivative of the matrix model potential W'(z) coincide to the order 
we have calculated, i.e. e\ = + 0(A 2N+4 ). This result was used in sec. 3.2 to compare the 
relation between a, and obtained from the matrix model with the one derived in ref. [31] 
using the M-theory curve (2.1c). The two results agree (see sec. 3.2). 



4.2.2 Method II: exact determination of the curve via Abel's theorem 

The condition imposed on the matrix-model curve by extremizing W e $ can also be discussed 
using the saddle-point solution, as in [7] and in sec. 4.1.2. 

The Riemann sheet structure of the generic cubic curve (4.30) is similar to that for 
the U(N)xU(N) model discussed in subsection 4.1.2, except that we are interested in the 
quotient of this curve by the involution z — + —z. 

In an eigenvalue basis the leading (sphere) part of the free-energy can be written 

F s = -Jd\ p(X)W(X) + JdX dX' [p(X)p(X') log(A - A') - \ P {X)p{X') log(A + A')] (4.45) 
and the sub leading (RP 2 ) part is 

Fry> = -\P J dAp(A)log(A) (4.46) 

where 

P(A) =9sJ2 5 ( X ~ ^) . «*>(*) = I dA ^T • ( 4 - 47 ) 

i Z — A 

We define 20 

s <=-hb {z)iz < (4 - 48) 

where Ai denote contours around the branch cuts near on sheet one. As before one can 
show that Si is the integrated density of eigenvalues along the cut near (so (4.48) is 

20 As in the U(iV) xXJ(N) theory, there are in general additional S"s corresponding to the other (possible) 
cuts of the curve. In particular, there is a variable So corresponding to the (possible) cut between sheets 
one and three around the point z = 

So = -— & u(z)Az 
27r * Jo 1 
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consistent with the perturbative definition Si = g s Mi) and 5p(X) = 5Si5(X — e.i). As in 
sec. 4.1.2 (up to terms which will not affect our discussion) 



dF s r 7 ! r 1 * , r 7 ! 

1 u — u — = / u 



dSi Ji 2 Ji 3 Ji 3 

jp p f 1 ,p f 1 p r 



u 

1 03 



(4.49) 



where 1$ denotes infinity on the ith sheet, and 0« is the point z = on sheet %. 
Next, we wish to extremize the effective superpotential (3.42) 

N dF 

= ~ E asf - 4F rp (4.50) 

(setting Ni = 1 and N = 0). By changing basis as in sec. 4.1.2 and varying w.r.t the 
coefficients of the arbitrary polynomials ri(z) and ti(z) in the matrix model curve (4.30), 
one obtains 

— (N — 2/3) [ h ( k -(N- 2(3) j h Ck + 2/3 ( k - 2/3 ( k = (modulo period lattice) . 

(4.51) 

This condition seemingly implies, by Abel's theorem, the existence of a function with a pole 
of order iV — 2(3 at I±, a zero of order N — 2(3 at J 3 , a pole of order 2(3 at 0i, and a zero of 
order 2(3 at O3, and regular everywhere else. However, there is one important caveat. In the 
undeformed {r\{z) = t\(z) = 0) curve (4.30), z = is a singular double-point. If a cut opens 
up between sheets one and three when r\(z) and ti(z) are turned on, then the points 0i and 
3 will be identical and the last two terms in (4.51) will not contribute, and the function will 
be regular at z = on all the sheets. 

(i) u(at) + m {fl = 1) 

We will now show that this problem has a solution: that for some choice of r±(z) and 
t±(z) in the matrix-model curve (4.30), there exists a function f(z) on this curve with divisor 
implied by (4.51) via Abel's theorem. For simplicity, we restrict ourselves to the special case 
in which N is even and o\ = 0; these conditions are such as to make W'(z) — W'(—z) a 
polynomial of order iV — 3. Consider r±(z) = and ti(z) oc z 2 (W'(z) — W'(— z)) 2 , i.e. the 
matrix-model curve 21 

(u - W'{zj) u{u- W'(-z)) = const x z 2 (W'(z) - W'(-z)) 2 . (4.52) 

The function 

defined on this curve, has the following asymptotic behavior near Jj, infinity on the three 
sheets: 

h ■■ f(z) ~ z N ~ 2 , h : f(z) ~ const , I 3 : f(z) ~ z~ N+2 (4.54) 



21 The polynomial ti(z) is of degree 2N — 2 or less only when only when N is even and o\ = 0. In the 
more general case, a more complicated choice of r\{z) and t\{z) would be necessary. 



26 



and the following behavior near Oj, z = on the three sheets: 



Oi : f(z) ~ z~ 2 , 2 : f(z) ~ const , 3 : /(z) ~ ^ (4.55) 

As in sec. 4.1.2 one can show that f(z) is regular everywhere else, and so satisfies the 
conditions implied by (4.51) and Abel's theorem. Assuming that the solution is unique 
(which we have not been able to prove), we find that extremization of W e s (via Abel's 
theorem) implies a matrix-model curve of the form (4.52). Upon redefining W'(z) oc P(z) — 
3A' N ~ 2 z 2 , this curve precisely agrees with the M-theory curve (2.7b), and the function (4.53) 
is proportional, up to a factor of z 2 , to y (2.6b). 

(ii) U(N) + Q(P = -1) 

To show that the problem has a solution, consider ri(z) = and t 1 (z) oc z~ 2 (W'(z) — 
W'(— z)) 2 , i.e. the matrix-model curve 

(u - W\z)) u(u- W'(-z)) = const x 1 ' J { ' (4.56) 

For this curve, there is a cut opening up at z = between sheets 1 and 3, so that the last 
two terms in (4.51) do not contribute. The function 

on the curve (4.56) has the following asymptotic behavior near Jj, infinity on the three sheets: 
h ■■ f(z) ~ z N+2 , h ■■ f(z) ~ const , h ■ f(z) ~ z~ N - 2 (4.58) 

and is regular near z = on all three sheets. As in sec. 4.1.2 one can show that f(z) is 
also regular everywhere else. Thus, the function (4.57) on the Riemann surface (4.56) has 
precisely the divisor specified by (4.51). Assuming that the solution is unique (which we 
have not shown), we find that Abel's theorem implies a Seiberg-Witten curve of the form 
(4.56). Setting W'(z) oc P(z), eq. (4.56) precisely agrees with the M-theory curve (2.7c), 
and the function (4.57) is proportional, up to a factor of z 2 , to y (2.6c). 



5 Seiberg-Witten differential from the matrix model 

In this section we will discuss the derivation of the Seiberg-Witten differentials for the Af = 2 
gauge theories described in section 2 from the matrix model point of view. 

As is the case for the S W curve, the SW differential can be determined order-by-order in 
A by using perturbative matrix model calculations. We will illustrate this in the first part 
of this section, after which we will very briefly discuss some other approaches, such as the 
one pursued in (version 3 of) ref. [7]. 
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5.1 Method I: perturbative determination of \$w 

On the first sheet of the Riemann surface, we have (for all the models) [4, 21, 3, 13] 

1 



X sw = zT(z)dz, T(z) = (tr 



. ; =E*- B - 1 <trtf B > 

> z 9)1 n =0 



(5.1) 



where (tr0 n ) is the gauge-theory vev of the adjoint field and T(z) is sometimes called h(z). 
The relation between gauge-theory and matrix-model vevs can be obtained using the methods 
in ref. [5] 



N 



M"> = £ 



i=l 



_d_ _d_ 

dS~ + dS, 



g s (tr$ n ) s2 



A? 



d 



(tr 0") = J2go.9s (tr $ n > 52 + 4(tr $") R p 2 

1=1 % 



(5.2a) 



(5.2b,c) 



where the first equation is for model (a), U(iV) xU(iV) with a bifundamental hypermultiplet, 
and the second equation is for models (b) or (c), U(iV) with a symmetric or antisymmetric 
hypermultiplet, respectively. It is understood that the rhs is evaluated at Si = (Si) (as well 
as Si = (Si) for model (a) and So = (So) for model (c)). Since we are only interested in the 
first two orders in perturbation theory, we may write [5] 



(tr*") = £ [M^ + ner'it^u) + |n(n-l) e r 2 (tr^> + 
i=i 

For model (c) we also have the extra terms 22 

M 5 nj0 + (tr^ 2 )5„ i2 + ... 



(5.3) 



(5.4) 



It follows from these expressions that for all three models, the leading term in the perturba- 
tive expansion for T(z) is given by 



N 1 

T(z) peit = E 



z - e, : 



(5.5) 



The matrix model expectation values (tr ^/a)s^ and (tr ^u)up 2 are given in eqs. (3.15), 
(3.48), (3.49) of sec. 3, and for all three models, we have 



{t^l)s* = 9 -j:% <tr* 2 )lRIp2 = 0. 



(5.6) 



In addition, for model (c), we will need 23 



(tr^ooW 2 



_^So_ 

2a R 



(5.7) 



22 Note that, for this model, (tr'J'oo) = because of the Sp-condition on <& o- 

23 Note that (trVE^)^ oc Sq does not contribute to (5.2b,c) because of the absence of a derivative w.r.t. So- 
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Using these results together with the above formula? leads to the following expressions for 
the one-instanton corrections (see section 3 for details about the notation) 



N 



•inst 



i l (z - ei) 2 
2 

+ 



2 T 1 T 

— ^ + E 

-^i j^i Rj e ij Ri j Rjhij 

Ti 



T 1 

~p2 I - 

-<4 j 'Hj 



(z - e,)3 it? 



(5.8a) 



inst 



2G ?: 



2 Gj 1 Gj G{ 1 

-Rj Rj^ij Ri j Rj9ij R% j 9ij R%&i 



+ 



(z - e,)3 i?2 



(5.8b) 



N 



inst 



2(-l) iV 



2 Gj 1 ^ Gj Gi 1 ^ 2Gj 

jyi Rj e ij Ri j Rj9ij Ri j 9ij Ri e i 



+ 



E 



+ 



2 Gj A_ 1 



rW 



(5.8c) 



where eqs. (5.8a), (5.8b), and (5.8c) correspond to models (a), (b), and (c), respectively. 
Next we use the identities (3.18), (3.53), and (3.58), together with the definitions (4.14), 
(4.39) to obtain 

T(z)dz = dlogP(^) - A N - 2 ^dK(z) + 0(A 2N - A/3 ) 



with 



P = 0, 
= 1, 

iN 



K(z) 



P(z) 2 



K{z) = (-1) 



N 



K(z) 



N 



z 2 P(-z) 
P(z) 2 

1 P(-z) 



3 1 



z 2 P(z) 2 z 2 P(z) 



(5.9) 

(5.10a) 
(5.10b) 
(5.10c) 



recalling that P(z) = Ui =1 (z - e[) and P(z) = Uti( z ~ 

These results are consistent with, using (5.1), what one obtains by expanding the M- 
theory result (2.2) on the first sheet of the Riemann surface given by (2.6a)-(2.6c). 



5.2 Other methods 

For the U(iV) xlJ(iV) gauge theory with a bifundamental hypermultiplet, we may combine 
eqs. (5.1), (5.2a), and (4.1) to derive the relation (on sheet one of the Riemann surface) 24 

N 

\ sw = zT(z)dz, T(z)=J2 
i=i 

24 In the equations in this section, it is understood that, after taking the derivatives, the results are to be 
evaluated at St = (Si) and Sj = (Si). 



_d_ _d_ 



u(z) 



(5.11) 
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where u(z) is the leading (sphere) part of the resolvent. Similarly, for the U(iV) gauge theory 
with one symmetric or antisymmetric hypermultiplet, we may combine eqs. (5.1), (5.2b,c), 
and (4.27) to derive (on sheet one) 



N Q 

X sw = zT(z)dz, T(z) = ^2ww-Uo{z) +Wi/i(z) 



(5.12) 



i=i 



where ojq{z) is the leading (sphere) part of the resolvent and ^1/2(2) is the subleading (RP 2 ) 
part. These results appeared (using a different approach) in [1, 19]), and earlier for the case 
of U(iV) without [50] and with [7] fundamental hypermultiplets. Now we may use eq. (4.3) 
to show, for U(iV)xU(A0, 



Since both Xsw and u are defined on all the sheets, this equation extends to the entire 
Riemann surface. A similar equation may be derived for U(iV) with a symmetric or anti- 
symmetric hypermultiplet starting from (5.12). 

Two methods for computing Xsw (hi addition to the approach used in sec. 5.1) present 
themselves. First, one may use perturbation theory to calculate the curve polynomials r\(z) 
and t\ (z) as functions of the S"s, and then use the curve equation to calculate the derivatives 
of u in eq. (5.13). In appendix A, we use the analog of this approach to determine the one- 
instanton contribution to T(z) for the U(iV) model with Nf < N fundamentals. Obtaining 
the one-instanton contribution for the theories considered in this paper is straightforward, 
but somewhat cumbersome, so we will not pursue it here. 

A second approach is to investigate the integrals around the A-cycles together with the 
behavior at infinity, and try to use this information to pin down the function T(z). This 
method was used in (version 3 of) [7] to determine T(z) for the U(iV) model with Nf < N 
fundamentals. For that case T(z) was given by where ip(z) was the function arising from 
Abel's theorem. Such a relation, involving the function f(z) given in (4.26), probably also 
holds for the U(N) xU(iV) theory. Recalling that f(z) (4.26) is proportional to y (2.6a), 
using (4.13), this implies T(z) = which is consistent with the M-theory result (2.2). For 
the U(iV) models the situation is less clear since f(z) (4.53), (4.57) and y (2.6b), (2.6c) differ 
by rescalings with z 2 . 

6 Summary 

This paper is part of a larger program aimed at studying the applicability of the Dijkgraaf- 
Vafa matrix model approach to M = 2 supersymmetric gauge theory theories, i.e. those 
theories amenable to the methods of Seiberg-Witten theory. Previously it has been shown 
that one can recover the ingredients of the Seiberg-Witten solution from the matrix model for 
theories with hyperelliptic Seiberg-Witten curves. In this paper we focused on three different 
models: (a) U(iV) xU(iV) with matter in the bifundamental representation, (b) U(iV) with 
matter in the symmetric representation, and (c) U(iV) with matter in the antisymmetric 




(5.13) 
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representation. Each of these theories is described by a cubic non-hyperelliptic Seiberg- 
Witten curve. Our goal was to determine the Seiberg-Witten curve and differential, as well 
as the order parameters and the prepotential, for the above models, entirely within the 
context of the matrix model, without reference to string/M-theory. Our results confirm the 
results previously obtained using M-theory. 

For models (a) and (b), a straightforward generalization of our earlier work (including a 
refinement of the prescription for Tjj, as discussed in [1]) produced expressions in complete 
agreement with earlier results in the literature. For model (c), U(iV) with matter in the 
antisymmetric representation, the naive extension of earlier work leads to discrepancies with 
previous results in the literature. The discrepancies occur when one expands around the 
simplest matrix model vacuum with IliU(Mj) gauge symmetry, where each Mi — > oo with 
Si = g s Mi fixed. However there are other vacua; one of these has Sp(Mo)x IliU(Mj) gauge 
symmetry [25], and hence an additional parameter So = g s M . When one expands around 
this vacuum, extremization of W eS (with N { = 1 and N = 0) leads to non-zero vevs for 
both Si and So, and all discrepancies are removed. What is missing is a guiding principle 
for the enlargement of the matrix-model vacua for the Q theory. (A similar enlargement of 
the vacua of the □□ theory seemingly leads to an inconsistency.) 

We close with an empirical observation about the form of the one-instanton contribution 
to the prepotential. For pure U(iV) gauge theory [5], or U(iV) with matter in the □ [7] or 
m (sec. 3.2.2) representations, one may verify, using our results, that the following equation 
holds 

2^-« = E^ (6-14) 

where the sum is over all the extrema of the superpotential W(z). Observe that the ex- 
pression (6.14) has a finite limit when the coefficient a multiplying the superpotential is 
taken to zero to restore M = 2 supersymmetry. For U(iV) xU(iV) with a bifundamental 
hypermultiplet (sec. 3.1.2), the sum also extends over the extrema of W(z): 

where the relative minus sign is due to the fact that W(z) enters with a minus sign in the 
superpotential (3.1). Finally, for U(iV) with Q (sec. 3.2.2), there is an additional contribution 
from the vacuum state at z = 0: 

^—^igj + ig, (««) 

However, at present we do not have an understanding of why these results are true. 
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Appendices 



A U(iV) + NfU curve and T(z) from perturbation theory 

In ref. [7], we derived the SW curve for the Af = 2 U(iV) gauge theory with Nf fundamental 
multiplets from the associated matrix model. Specifically, the saddle-point solution of the 
matrix model gives rise to the equation for the resolvent 25 

u 2 {z)-W'{z)u{z) + \f(z)=0 (A.l) 

where W'(z) = a~[\f =1 (z — e*) and using (4.7) 

/(*) = 4g.(tr ~ ^'^ ) ) = 4 E cn(z) g.(trm) (A.2) 

is an (N — l)th order polynomial. Defining y(z) = — 2u(z) + W'(z) one may rewrite eq. (A.l) 
as a hyperelliptic curve 

y 2 = W'(z) 2 -f(z). (A.3) 

In ref. [7], we determined the polynomial f(z) using Abel's theorem (see also [4]). In 
this appendix, we will show how this polynomial can be evaluated using the perturbative 
solution. The lowest order contribution to the matrix model vev is [5, 7] 

g s (tr($ n )) = aA 2N ~ N f ]T ^e™ + 0(A m ~ 2N f) (A.4) 

where = llj=i(ei + mi) and Ri = Yl^fe — ej). It follows from (A.4) that 

f(z) = 4aA 2N ~ N f [W'(z) - W'ia)] £ D , Li v + 0(A* N ~ 2N f) . (A.5) 

; Ri{z-ei) 

Now using the fact that W'(ei) = together with 

v Li _ Ilfl 1 (z + m I ) f 

V nf^-e.) T(2) (A - 6) 

where T(^) is the polynomial part of Ylfl^z + mi)j Yu=\k z ~ e i-> we finally obtain 

/(z)=4a 2 A 2 ^/ ( IKz + m^-f^lKz-eM + 0(A 4N - 2N f) (A.7) 

which is precisely the result obtained by Abel's theorem in ref. [7]. 

When Nf < N, the polynomial T(z) vanishes. Moreover, in that case, eq. (A.7) is exact 
to all orders in A [7], which points to the existence of a non-renormalization theorem. 



25 



We have converted to the notation in this paper; to convert back let u(z) — > — Su(z). 
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Similarly, in (version 3 of) [7] (see also [50]), we derived the expression for the Seiberg- 
Witten differential for this theory (with Nf < N) from the saddle-point approach 



N dv 

X sw = z T{z) dz, T{z) = -\ E ^ 

1=1 1 



(S) 



+ 



y - W'(z) f{z) 
2y f(z) ' 



(A.8) 



Let us now evaluate this expression using only the perturbative solution. As in eq. (4.8), we 
have 



f(z)=4W'(z)J2-^- + 0(Sf) 

Z di 



(A.9) 



so that, using eq. (A. 3) 



N 



dy 



N 



df 



N 



W'(z) 



+ O(Si) 



W"(z) 

y 



+ 0(Si) . 



(A.10) 



Evaluating this expression at Si = (Si), the O(Si) term is subleading in A, so the leading 
order contribution to T(z) is 



T(z) 



W"{z) + \{y-W\z)) 



(A.11) 



However, as we showed in [7], this expression is exact, again pointing to the existence of a 
non-renormalization theorem (when Nf < N). 

B Some technical details 

In this appendix we collect some technical details of the perturbative calculations performed 
in section 3. 



B.l V(N)xU(N) 

The integration over the quadratic fields of the matrix-model partition function (3.3) yields 
(up to an ej-mdependent quadratic monomial in the Si, S^s) 



N 



N 



F S (S,S) 



Si 



- g + hY. s! >og {^) + 1 E s,s, io g 

% _L % A. % A. -f—% 

N N 



aRiA 2 , 
Si 



+ 



-aRiA 2 / i=irfi 

EE bg + 

i j \ A / rt>3 



N 

+ EE^iog( : f ) + 



(B.l) 



Cj, hij = ei — Cj and Ri = flj^e^. The term F^ n \S, S) is an nth order 
polynomial in Si and Si arising from planar loop diagrams built from the interaction vertices. 



where = e 
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The contribution to F S (S, S) cubic in Si and Si, 

2 



aF&(S,S) = 



S: 



1 



Sf 



3 E R . ( E .. ) 4E R .E^ 



+ 2 EEE 



SiSkSe 

■ RiGik&i£ 



k^i i^i,k e ik e it 



i kjti - rt « e «fc tjki tit 
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/ 1 , i , , . j\i I Oik n-^i e i£ i ' '• i • 1 , 1 



sfs k 



i k^i AL ^ik i k 



i k ± ^' u ik 



„ \ ^ \ -> \ -« SiSkSe ^ x -> x x SiSkSi 



, (B.2) 



is obtained by adding to the result from ref. [5] the planar two-loop diagrams containing B 
loops that can be drawn on a sphere. 

Using the above expressions to calculate (3.10) and extremizing (3.11) leads to 



(Si) 



Ri 



A 



271; 



+ 



2T k 



271 



+ 



27/ 



+ EE 



37^. 
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key 
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+ EE 
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+ E 
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RkRl&klhil) 

Ti Tk 



+ EE 

k i 

fk 



T k 



RiRkhikhu R\hikhik RkRehieh 



k£ , 



+ 0(A 2N ) 



Rfhik Rlhfk RiRkhj 



ik / 



(Si) 



Ri 

rN 



1 + A J V above expression with <-> > + 0(A 



2N\ 



(B.3) 



where Tj = n,=i( e i — <%•), Tj = fT - =1 (ei — e^), Ri = Wj^e-i — ej), and various constants have 
been absorbed into a redefinition of the cut-off A. 

The gauge-coupling matrix (3.20) can be expanded as 



Tij T^j 



pert L j^Nd^d) 
d=l 



(B.4) 



and similarly for T\~ 3 and Tjj. Using (B.l), (B.2), (B.3), and then reexpressing the result in 
terms of dj and hi using (3.16), the first two terms in the expansion (B.4) can be determined. 
The perturbative contribution is (up to additive constants) 



2nirff rt = 



■2£log 

k+i 



<H - CLk 
A 



Eiog 

k 
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2 log 
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(B.5) 



and the one-instanton contribution, after some algebra, is 



2-nirf 
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(B.6) 
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where now = a* - a j: Kj = a* - Ri = Uj=a a>ij, R% = Ylj^ii^i ~ T i = ITjii(aj - 
and Tj = n^life — %)■ The expressions for are obtained from those above by letting 
di <-> 5j. Finally, 

2^ ert = -bg(*^) (B.7) 



A 



(l) _ ^ / 2 y-v 1 | 1 | | 2VV I (i < ) j- a < > a ^ 



These expressions may be expressed succinctly as the second derivative of the prepotential 
(3.22). 

B.2 \]{N) with m or g 

Using the Feynman rules, the cubic contribution to the sphere part of the free energy can 
be shown to be 

«f?\s,§) = f£f (zf) -|Ef EE-^EEf^Ef 

CCC C2 Q C2 C 1 

+ 2 z^ z^ z^ p ~ ~ z^ z^ d p 2 + z^ z^ p ~ z^ ~ 

j k^i i^i ^i^iktil i ^t^ik i k ^HiM ^il 
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i fc^i £ ^ikyit i k e n iyikyie i k n i9ik 

CCC CCC CC 2 C C 2 

2 E E yv OpOjOj ZZ ° °i 1 'JQ^t 

+i 2 ~ gjj + 5*5- terms + 5 3 terms (B.8) 

where e„ = e, t - ej, = e { + e,, i?, = Uj^i^ij, and i? = — e«) Ej(V e i)- Tnis result 
was obtained by adding to the result from ref. [5] the new planar two-loop diagrams that can 
be drawn on a sphere and contains Y and ghost loops. Similarly the quadratic contribution 
to the PJP 2 part of the free energy can be shown to be 



c 2 1 cc cc c 2 cc cc 

Z^ ~ Z^ r? ~ ~ ^ 2L 



R%&i dil j-^i Ri&i&ij j Ri&i9ij Ri^-% Ri&i R()&i 

+S'£ terms, (B.9) 

The terms in (B.8) and (B.9) involving Rq come from diagrams containing \l/oo legs. Since 
^00 £ sp(M ) for /3 = —1, one must use (1/2R ) {5^51 + J ac Jbd) for the propagator [51]. 
Observe that eqs. (B.8) and (B.9) obey the relation given in footnote 14. 

Using the above expressions to calculate (3.42) and then (3.43), the solution for (Si) can 
be evaluated in an expansion in A 
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(B.10) 



where Gi = ^fY^jipi + e j)- To evaluate (3.61) perturbatively as 



(B.11) 



d=i 



we use eqs. (3.40) and (B.8) in eq. (3.61), evaluate the resulting expression using eqs. (3.46) 
and (B.10), and then use the results of section 3.2.1 to re-express the entire expression in 
terms of a,. The perturbative contribution is (up to additive constants) 
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and the one-instanton contribution, after some algebra, is 
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(B.13) 



+ (i <-> j) 



where now = aj — a^, = + a^, i?j = Ylj^i a ij, an d Gj = Ylj( a i + Equations 
(B.12) and (B.13) may be written succinctly as the second derivative of the prepotential 
(3.62). 
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